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CHAPTER  1 

INTRODUCTION 


Background 


Continued  research  of  the  subsonic  boundary  layer  laminar-to-turbulent  transition 
zone  wall  pressure  fluctuations  has  resulted  in  the  formulation  of  an  empirical  space-time 
correlation  model  of  the  pressure  field.  Josserand  and  Lauchle  ^  developed  the  empirical 
model,  by  determining  the  parameters  of  assumed  functional  forms,  from  extensive 
measurements  of  spot  formation  in  the  transition  zone.  They  measured  a  "zero-one" 
function,  an  indicator  function,  that  identifies  the  presence  of  turbulent  spots.  The  space- 
time  correlation  model  of  the  indicator  function  provides  the  fundamental  statistics  of 
the  stationary,  non-homogenous  pressure  field.  However,  the  preferred  format  of  this 
pressure  field,  for  descriptive  and  interpretation  purposes,  is  the  wavevector-frequency 
spectrum.  The  wavevector-frequency  spectrum  is  obtained  by  Fourier  transformation  of 
the  space-time  correlation.  Two  complications  arise  when  attempting  the  Fourier 
transformation.  First,  the  relatively  complex  form  of  this  empirically  derived  correlation 
function  defies  the  analytical  integration  required  to  produce  the  desired  wavevector- 
frequency  spectrum.  Second,  since  the  pressure  field  is  non-homogeneous,  the  resulting 
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wavevector-frequency  spectrum  will  be  a  function  of  reference  position.  This 
dependence  on  reference  position  precludes  a  concise  description  of  the  field.  That  is, 
for  a  homogeneous  pressure  field,  a  single  wavevector-frequency  spectrum  describes  the 
field.  For  a  non-homogeneous  pressure  field  no  single  wavevector-frequency  spectrum 
describes  the  field.  In  fact,  an  infinity  of  spectra  are  required.  Therefore,  although  the 
wavevector-frequency  spectrum  is  the  preferred  format,  a  more  concise  spectral 
description  is  desired.  Strawderman  has  recently  identified  a  full  spectral  definition  of  a 
non-homogeneous  field.  This  thesis  will  investigate  the  application  of  both  spectral 
forms  to  the  study  of  boundary  layer  transitional  flow  pressure  fluctuations. 


Goals  and  Scope  of  Work 


The  objective  of  this  research  is  to  extend  the  work  of  Josserand  and  Lauchle  by 
computing  two  forms  of  the  wavevector-frequency  spectrum  of  transition  flow.  First,  the 
space-varying  wavevector-frequency  spectrum  will  be  estimated  using  computer 
implemented  numerical  methods.  Second,  in  an  attempt  to  provide  a  more  concise 
wavevector-frequency  description  of  this  non-homogeneous  field,  a  recently  developed 
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spectral  format  named  the  "Two  -  Wavevector  -  Frequency  "  spectrum  will  be 
computed.  These  computations  will  require  multi-dimensional  Fourier  transformation  of 
the  empirically  derived  space-time  correlation  function  of  the  boundary  layer  indicator 
function.  The  extent  of  the  fundamental  dimensionality  of  these  transformations  is  great. 
In  general,  one  temporal  and  two  independent  spatial  vectors  are  required.  Techniques  to 
image  these  dimensions  are  not  generally  available  and  are  beyond  the  scope  of  the 
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present  investigation.  Therefore  it  is  necessary  to  select  a  useful  but  reduced  domain  of 
interest. 

It  is  known  that  the  major  features  of  turbulent  flow  are  contained  in  the 
streamwise  component.  Investigation  of  the  characteristics  of  the  space-time  correlation 
of  transitional  flow  reveal  similar  major  features.  Therefore  it  is  considered  appropriate 
to  limit  the  dimensionality  of  the  current  investigation  to  the  streamwise  component  of 
transition  flow. 


This  will  require  satisfying  the  following  goals; 

1.  Implement  a  computer  model  of  the  indicator  function  space-time  correlation 
function. 

2.  Compute  the  Fourier  transform  of  the  indicator  function  space-time  correlation 
function  on  streamwise  wavenumber  and  frequency,  at  selected  reference 
positions  in  the  transition  zone,  to  estimate  the  space  varying  wavenumber- 
frequency  spectrum  of  transitional  flow. 

3.  Compute  the  Fourier  transform  of  the  indicator  function  space-time  correlation 
function  on  streamwise  wavenumber,  frequency,  and  reference  position  to 
estimate  the  two  wavenumber-frequency  spectrum  of  transitional  flow. 
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CHAPTER  2 


SPACE- VARYING  WAVENUMBER-FREQUENCY  SPECTRUM  OF 
TRANSITION  FLOW  PRESSURE  FIELD 


The  wavenumber-frequency  spectrum  of  boundary  layer  transitional  flow  wall 
pressure  fluctuations  is  defined  ^  as 

j  +—  +— 

(27?)  _  (1) 

Here  r|^  denotes  absolute  position  within  the  transition  zone,  e,  is  longitudinal  spatial 
lag,  relative  to  77,,  and  £3  is  transverse  spatial  lag.  Temporal  lag  is  denoted  by  r. 
Longitudinal  and  transverse  wavenumber  are  denoted  by  and  /.'j,  co  is  radian  temporal 
frequency. 

The  function  R,[r|^,£^,£J,r)  is  the  space-time  correlation  of  the  indicator 
process  and  /?„^(£,,£3,t)  is  the  space-time  correlation  of  the  turbulent  flow  wall 

pressure  fluctuations,  presumed  homogeneous  within  the  turbulent  spots  of  the  transition 
zone.  The  correlation  function./?/(r7,,£, ,£3, t)  is  derived  empirically  from  measurements 

of  the  statistics  of  transition  zone  spot  formation  and  contains  the  non-homogeneous 
characteristics  of  interest  here.  The  correlation  ^;„,i(£,,£3,T)  is  the  homogeneous 

component  and  contains  no  information  from  the  transition  zone.  Considering  this,  and 
recalling  the  previous  discussion  regarding  the  desire  to  reduce  the  dimensionality  to 


streamwise  flow  only  (Cj  =  0 ),  it  is  appropriate  to  write  the  pertinent  wavenumber- 
frequency  spectrum  as 
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(2) 


Space-Time  Indicator  Correlation  Function 


O 

The  space-time  correlation  of  intermittency  has  been  determined  empirically 


to  be 
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Here,  7u  is  the  intermittency  measured  by  the  upstream  probe  and  is  the 
intermittency  measured  by  the  downstream  probe  .  The  convection  velocity  is  Ax 
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is  the  streamwise  length  of  the  transition  zone,  A,  is  the  normalized  longitudinal 
separation  inside  the  transition  zone,  e,/ Ax,  and  A  =  -ln[l-exp(-4.27/ Aj:)]  with  Ax  in 
meters.  In  general,  the  intermittency  is  written  as  y{z)  =  1  -  exp[-(l  -  3.4z)z^]  with 

z=  Tji/Ax. 

It  is  useful  to  identify  the  t]^  dependence  of  (3)  for  £3  =  0.  This  dependence 

is  qualitatively  demonstrated  in  the  contour  plots  shown  in  Figures  1, 2,  and  3.  These 
contour  plots  were  generated  for  Ax  =  1.04  m  and  three  arbitrarily  selected  values  of  77,: 
0.1040,  0.5200,  and  0.8320  meter  from  which  z  ranges  from  0.1  to  0.8.  Note  that  each  of 
the  figures  can  be  described  in  terms  of  only  two  processes  that  occur  to 
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TT  seconds 


Figure  1.  Contour  Plot  of  Space-Time  Correlation  of  Intermittency;  z  -  0.8 
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Figure  2.  Contour  Plot  of  Space-Time  Correlation  of  Intermittency;  z  =  0.5 

varying  degrees  as  a  function  of  reference  location  ( 77,).  One  process  has  correlation  that 
varies  only  with  spatial  lag.  The  other  process  is  highly  correlated  over  only  a  small 
temporal  region.  The  time  lag,  associated  with  the  peak  in  correlation,  varies  with  spatial 
separation  indicating  propagation.  This  propagating  component  has  a  velocity  of 
approximately  8.23  m/sec,  and  is  consistent  with  the  selected  convection  velocity. 

In  summary,  these  figures  demonstrate  the  expected  convective  component  of 
intermittent  turbulent  flow,  as  well  as  the  position  dependent  (non-homogeneous)  nature 
of  the  transition  flow  pressure  field,  in  terms  of  space-time  correlation. 
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Figure  3.  Contour  Plot  of  Space-Time  Correlation  of  Intermittency;  z  =  0.1 


Finite  Fourier  Tran.sform  for  Computer  Estimation  of 
Tran.sition  Flow  Wavenumber  vSpectrum 


The  pertinent  wavenumber  spectrum  of  transition  flow  is  defined  by  (2). 

However,  as  previously  staled,  a  digital  computer  will  be  used  to  generate  numerical 
spectral  estimates.  The  computer  implementation  requires  a  finite,  discrete,  form  of  (2). 
The  consequence  of  the  required,  finite  length,  space-time  sampling,  is  spectral  amplitude 


bias,  where  spectral  amplitude  bias  is  defined  ^  as  the  difference  between  the  spectral 
estimate  and  the  true  spectrum. 

The  concept  of  spectral  bias  can  be  developed  from  the  realization  that  finite, 
discrete  space-time  correlation  sample  sets  are  obtained  from  continuous,  infinite  space- 
time  correlation  functions,  by  multiplication  by  appropriate  sampling  functions.  The 
wavenumber-frequency  spectrum  is  therefore  the  Fourier  transform  of  a  product.  The 
Fourier  transform  of  a  product  is  equivalent  to  the  convolution  of  the  Fourier  transforms 
of  the  individual  functions.  Therefore,  the  streamwise  wavenumber-frequency  spectrum 
may  be  thought  of  as  the  convolution  of  the  two  dimensional  Fourier  transform  of  the 
space-time  sampling  function  with  the  true  wavenumber-frequency  spectrum.  The 
Fourier  transform  of  a  finite  length  space-time  sampling  function  will  contain  a  main 
response  lobe  of  finite  width.  Therefore,  a  spectral  estimate  at  any  given  wavenumber 
and  frequency,  comprises  an  integral  of  the  true  wavenumber-frequency  spectrum  giving 
rise  to  spectral  amplitude  bias. 

Spectral  bias  is  the  combined  effect  of;  1)  aliasing^  resulting  from  sampling  at 
uniform  intervals,  2)  decreased  spectral  resolution^  resulting  from  finite  length  sample 
sets,  3)  sidelobe  leakage"^  caused  by  discontinuities  at  the  sample  set  boundary. 
Windowing  techniques  exist,  and  will  be  used  in  this  research,  to  mitigate  the  effects  of 
sidelobe  leakage.  Spectral  aliasing  is  an  inherent  consequence  of  sampling.  However, 
aliasing  occurs  in  a  systematic  manner  allowing  meaningful  data  interpretations  to  be 
made  for  some  relatively  uncomplicated  cases. 

Decreased  spectral  resolution  is  particularly  troublesome.  The  spectral 
amplitude  bias  caused  by  decreased  spectral  resolution  can  be  reduced  or  eliminated  by 
appropriate  scaling  terms- if^the  resolution  and  bandwidth  of  the  process  of  interest  are 
known  apriori.  The  resolution  can  be  determined  from  knowledge  of  the  space-time 
sampling  parameters.  In  as  much  as  a  goal  of  this  research  is  to  provide  unprecedented 
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images  of  two  transitional  flow  pressure  field  spectral  formats,  it  follows  that  the  process 
bandwidth  is  a  desired  result,  and  therefore,  not  known  apriori.  Consequently,  any 
computer  implemented  spectral  calculation  is  subject  to  spectral  amplitude  bias  and  thus 
can  only  estimate  the  true  spectrum^.  It  should  be  emphasized  that  the  impact  of  spectral 
bias  is  to  introduce  an  unknown  amplitude  scaling  term.  Spectral  shape  characteristics, 
which  provide  the  fundamental  basis  for  interpretation,  are  not  affected.  Therefore,  the 
impact  of  spectral  bias  on  this  research  is  minimal.  Considering  these  factors  it  is 
deemed  appropriate  to  write  the  required  finite,  discrete  form  of  (2)  as 

^  N  M 

'F/(r7,,^,,(y)«  (/iAe,)iv^(mAT)i?;(77,,nAe,,mAT) 

/.=!  m  =  l  (4) 

where  the  proportionality  is  used  to  reflect  the  amplitude  uncertainty. 

This  implementation  requires  selection  of  finite  sets  of  discrete  values  for  the 
transform  variables  n,  m,  Ae,  and  At  as  well  as  spatial  and  temporal  weighting 
functions  w^  (nA£,),and  w^(mAt).  Since  /?y(T7, ,£,, £3, t)  was  empirically  derived  it  is 

considered  appropriate  to  select  the  required  computer  model  sampling  parameters 
consistent  with  those  of  the  measurement  apparatus.  Pertinent  parameters  of  that 
measurement  apparatus  are  discussed  in  the  following  section. 

The  space-time  correlation  of  yansition  flow  pressure  was  obtained  by  indirect 
methods.  Measurements  of  boundary  layer  spot  formation  in  the  transition  zone  were 
obtained  from  an  array  of  hot  film  sensors^  as  shown  in  Figure  4.  The  actual 
measurement  process  consisted  of  sequential  data  acquisitions  from  discrete  sensor  pairs. 
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Figure  4.  Probe  location  array 

This  data  set  was  used  to  generate  parameters  for  assumed  functional  forms  of  the 

indicator  function  space-time  correlation.  However,  noting  that  the  data  were  obtained 
from  an  "array"  of  spatially  distributed  sensors,  it  is  useful  to  perceive  of  R,[ri^,£^,£^,x) 

having  been  obtained  by  direct  measurement.  That  is,  consider  replacing  the  hot  film 
sensors  by  pressure  sensors.  Acquiring  data  from  all  pressure  sensors  ,  simultaneously, 
would  provide  a  direct  measurement  of  the  transition  zone  wall  pressure  fluctuations. 
Array  signal  processing  techniques  could  then  be  applied  to  the  discrete  pressure 
measurements  to  produce  the  wavenumber  spectrum  of  transition  zone  wall  pressure. 

This  concept  of  considering  the  hot  film  sensors  to  be  a  "virtual  array"  of  pressure  sensors 
is  the  basis  for  the  selection  of  the  sampling  characteristics  of  the  computer  model. 
Because  only  the  streamwise  component  is  considered  in  this  current  research,  the 
required  transform  parameters  are  the  rate  and  quantity  of  temporal  and  spatial  samples 
from  the  elements  of  a  linear  array. 

The  temporal  sampling  rate  was  selected  to  preclude  temporal  aliasing  at 
frequencies  of  interest  consistent  with  those  in  reference  (1).  This  requires^®  at  least  two 
samples  of  each  period  at  the  highest  frequency  of  interest.  Typically,  frequencies  of 
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100  Hz  or  less  were  considered  requiring  a  Nyquist,  or  minimum,  sampling  frequency  of 
200  Hz.  Sampling  in  excess  of  this  minimum  value  is  desired.  Therefore,  the  temporal 
sampling  rate  was  selected  to  be  250  Hz  corresponding  to  a  sample  interval  of  0.0040 
seconds.  Associated  with  this  temporal  sampling  rate  is  the  assumption  that  the  measured 
data  was  low  pass  filtered  to  eliminate  frequencies  greater  than  100  Hz.  The  number  of 
temporal  samples  was  selected  to  be  513. 

Similar  sampling  criteria  applies  to  the  spatial  domain.  Here,  the  desire  to 
simulate  the  measurement  apparatus  defines  the  spatial  sampling  rate.  From  Figure  4,  the 
linear  section  consists  of  thirty  two  locations  separated  by  2.54  cm,  thus  dictating  a  2-54 
cm  spatial  sampling  rate.  Strawderman^  1  .identifies  the  consequence  of  this  uniform 
spatial  sampling  rate  as  spatial  aliasing  of  spectra  at  wavenumbers  greater  than  ir/Ae,. 
This  means  that,  for  the  2.54  cm  spatial  sampling  rate,  components  of  the  wavenumber 
spectrum  in  the  range  !/:,!<  123.7  rad/m  will  be  located  at  their  actual  (unaliased) 
position.  However,  spectral  components  having  higher  wavenumbers  will  be  aliased  into 
the  wavenumber  range  of  I  k^  k  123.7  rad/m.  The  true  wavenumber  of  unaliased  spectral 
components  can  be  estimated  by  realizing  that  the  true  wavenumber  location  is  separated 
from  the  aliased  wavenumber  by  an  integer  multiple  of  2  %/  A£^  Therefore,  although 

the  occurrence  of  spatial  aliasing  associated  with  implementing  uniform  spatial  sampling 
confounds  data  interpretation,  it  does  not  preclude  it. 

The  nominal  length  of  the  measurement  array  is  approximately  0.81  m. 
However,  from  Table  1,  the  length  of  the  transition  zone  was  always  somewhat  greater 
than  the  array  length  (1.00  to  1.24  m).  Therefore,  in  an  attempt  to  simulate  the  spectral 
resolution  of  an  aperture  consistent  with  the  nominal  length  of  the  transition  zone,  but  not 
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deviate  greatly  from  that  of  the  actual  measurement  apparatus,  the  number  of  sensors  was 
selected  to  be  forty  one. 


Table  1.  Parameters  of  the  Measured  Transition  Zones 


1 


Flow  Velocity 

Measured 

Transition  Length 

Transition  Point 

t/_  (m/s) 

Intermittency  Range 

Ax  (m) 

Xo  (m) 

11.2 

0.08  -  0.60 

1.24 

0.990 

11.77 

0.17-0.80 

1.04 

0.914 

12.6 

0.40  -  0.99 

1.00 

0.863 

Implementation  of  the  spatial  and  temporal  sampling  parameters  implies 
multiplication  of  by  finite  length,  unit  amplitude  sampling  functions.  The 

Fourier  transform  of  these  finite  length,  unit  amplitude,  sampling  functions  are  infinite 
length,  periodic,  spectral  representations  that  comprise  mainlobe  and  sidelobe  structure. 
Unit  amplitude  weighting,  applied  to  a  given  length  sample  set,  produces  the  narrowest 
mainlobes  and  highest  sidelobes.  The  narrow  mainlobe  is  advantageous.  However,  the 
high  sidelobes  typically  contribute  unacceptable  spectral  estimation  error  via  sidelobe 
leakage.  This  sidelobe  leakage  can  be  reduced  by  the  use  of  window  functions  having 
other  than  unit  amplitude  coefficients.  Many  windows  have  been  developed.  In  general, 
the  consequence  of  other  than  unit  amplitude  windows  is  broader  mainlobes  and  reduced 
sidelobes.  Therefore,  the  decision  regarding  the  need  for  a  window  involves  a  trade-off 
between  spectral  resolution  and  sidelobe  leakage.  Because  this  research  involves  the 
pressure  field  of  a  turbulent  boundary  layer,  known  to  contain  high  energy  levels 
concentrated  at  specific  regions  of  the  wavenumber  domain,  a  window  was  deemed 
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necessary  to  reduce  sidelobe  leakage.  A  Taylor^^  weighting  function  was  selected,  for 
both  the  spatial  and  temporal  windows.  This  weighting  function  may  be  written  as: 


w{p)  = 
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/i-I 


F(0)  +  2^  F{m)  cos(mp) 
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and  n=  number  of  "near-in"  nulls,  =  , 


A  =  -\nl^R  +  4R^y 


V  2 


!£1 

/?  =  10“ ,  and  S=the  maximum  sidelobe  level  in  dB. 


(5) 

(6) 

(7) 

(8) 


Examples  of  Estimated  Space-varving  Transition  Flow 
Wall  Pressure  Wavenumber  Spectra 


The  space-varying  wavenumber  spectrum  was  calculated  for  two  frequencies 
at  reference  positions  of  z  =  0.1, 0.5,  and  0.8.  These  reference  positions  are  consistent 
with  those  chosen  for  the  space-time  correlation  of  intermittency.  The  computed  spectra 
for  50  Hz  are  shown  in  Figure  5.  Note  that  the  abscissa  is  wavenumber,  non- 
dimensionalized  by  transition  zone  length  (1.04  m  for  this  case).  Interpretation  of  these 
spectra  reveal  that  there  are  clearly  two,  and  perhaps  three,  coherent  processes  visible 
above  the  underlying  incoherent  component.  The  most  prominent  is  that  occurring  at 
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Figure  5.  Estimated  Space-Varying  Streamwise  Wavenumber  Spectrum;  50  Hz 

zero  wavenumber.  Energy  in  this  region  of  the  wavenumber  spectrum  can  occur  from  at 
least  two  different  processes.  One  possibility  is  that  the  spectra  are  of  acoustic  origin. 
Another  is  that  they  reflect  the  mean  values  of  the  random  processes.  The  origin  of  this 
component  was  resolved  by  eliminating  the  mean  value  from  the  space-time  correlation 
function  and  re-computing  the  space-varying  wavenumber  spectra.  As  shown  in 
Figure  6,  the  prominent  feature  at  ki=0  has  been  eliminated.  Clearly  this  component  was 
the  result  of  the  mean  value  contained  in  the  space-time  correlation  function.  In  as  much 
as  this  research  focuses  on  the  dynamic  characteristics  of  transition  flow  wall  pressure, 
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Figure  6.  Estimated  Space-Dependent  Streamwise  Wavenumber  Spectrum;  50  Hz 

and  the  fact  that  the  mean  value  could  preclude  observation  of  these  characteristics,  all 
following  estimates  of  wavenumber  spectra  presented  here  will  be  computed  from  zero 
mean  estimates  of  the  correlation  functions.  The  second  coherent  component  observed 
in  Figure  6,  occurs  at  a  non-dimensional  wavenumber  of  approximately  38.  This 
wavenumber  is  associated  with  a  phase  velocity  of  8.23  m/sec  suggesting  that  this  is.the 
convective  component  of  the  turbulent  spots.  This  component  is  most  visible  when  z  = 
0.1,  and  0.8.  At  z  =  0.5  it  is  substantially  masked  by  the  increased  amplitude  of  the 
incoherent  component.  Note  that  the  amplitude  of  this  convective  component  varies  with 
reference  position,  and  as  expected,  is  therefore  non-homogeneous. 
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The  third  coherent  component  is  only  weakly  visible,  at  a  non-dimensional 
wavenumber  of  approximately  88,  for  z  =  0.1  and  0.8.  This  process  is  completely 
masked  by  the  incoherent  process  at  z  =  0.5.  It  is  interesting  to  note  that  the  phase 
velocity  associated  with  this  wavenumber  ( 3.6  m/sec),  is  approximately  half  that  of  the 
convective  component.  The  importance  of  this  observation  will  be  developed  in  a  later 
discussion. 

The  incoherent  process  is  revealed  as  the  relatively  constant  amplitude  over 
the  wavenumber  domain.  Note  that  this  process  varies  with  z,  so  it  is  therefore  non- 
homogeneous  with  the  greatest  amplitude  occurring  at  z  =  0.5.  This  is  expected  because 
the  maximum  temporal  turbulent  spot  generation  rate  occurs  in  the  center  of  the  transition 
zonc^^. 

Figure  7  shows  the  estimated  space-varying  wavenumber  spectra  at  100  Hz. 
Interpretation  of  these  spectra  reveal  that  the  same  processes  that  exist  at  50  Hz,  exist  at 
100  Hz  and  are,  to  some  extent,  more  discernible.  The  second  convective  component 
remains  non-homogeneous  and  occurs  at  a  non-dimensional  wavenumber  of 
approximately  76,  consistent  with  the  increase  in  frequency.  A  third  component  is  seen 
at  a  wavenumber  of  approximately  -116.  This  wavenumber  is  likely  the  alias  of  the 
actual  wavenumber.  This  hypothesis  results  from  the  following  considerations.  First 
note  that  the  wavenumber  of  interest  (-116)  implies  that  the  direction  of  propagation  is 
opposite  that  of  the  previously  identified  convective  components  and  therefore  not 
consistent  with  any  expected  phenomena.  A  physically  possible  phenomena  is  spatial 
aliasing.  From  the  previous  discussion  we  know  that  spatial  aliasing  will  occur  at  integer 
multiples  oil  id  A£, .  The  estimate  of  the  true  wavenumber  is  determined  as 
-116  +  2  7C/0.0254  =  131.  This  wavenumber  suggests  a  phase  speed  of  4.8  m/s. 

Therefore,  this  wavenumber  produces  a  reasonable  interpretation  that  is  consistent  with 
expectations  in  both  wavenumber  magnitude  and  direction  of  propagation.  This  phase 
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Figure  7 .  Estimated  Space-Dependent  Streamwise  Wavenumber  Spectrum;  100  Hz 
speed  has  increased  somewhat  with  frequency  but  is  still  approximately  half  that  of  the 
convective  component.  Specifically,  from  Figure  6,  the  phase  speed  of  the  third 
component  is  0.3 1 1/_,  while  from  Figure  7  the  phase  speed  is  0.41  f/_. 

Lauchle^^  has  reviewed  the  results  of  several  studies  that  identify  the 
characteristics  of  turbulent  spot  development  in  a  transitioning  boundary  layer.  He  notes 
that  turbulent  spots  typically  demonstrate  a  growth  in  spatial  extent  with  streamwise 
propagation.  This  spatial  growth  rate  implies  that  the  celerity  of  the  leading  and  trailing 
evolving  turbulent  spot  edges  are  different  It  is  noted  that  Mautner  and  Van  Atta 
determine  leading  and  trailing  edge  pressure  convection  velocities  to  be  0.82ue  and 
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0.55ue  where  "ue  is  the  velocity  at  the  outer  edge  of  the  laminar  boundary  layer".  This 
edge  velocity  is  typically  0.99  for  a  zero-pressure  gradient  flat  plate  flow.  From 
Figure  7  a  similar  data  comparison  is  0.7  and  0.41  respectively.  Although  not 
exact*  ,  these  values  appear  consistent  with  the  reported  characteristics.  Therefore,  it 
seems  reasonable  to  conjecture  that  the  features  observed  in  the  space-varying 
wavenumber  spectra  are  the  result  of  the  turbulent  spot  generation  process.  The 
incoherent  process  is  seen  to  be  non-homogeneous. 


vSummarv  of  the  Space-varying  Wavenumber  Spectra. 


The  utility  of  the  space-varying  wavenumber  spectrum  is  apparent  from  only 
these  few  examples.  The  examples,  which  constitute  unprecedented  wavenumber- 
frequency  spectral  estimates  of  the  transitional  flow  pressure  field,  demonstrate  that 
spectral  decomposition,  in  the  wavenumber  domain,  allows  separation  and  identification 
of  processes  not  obviously  apparent  in  the  space-time  correlation  function.  The  examples 
also  demonstrate  the  position  dependent  nature  of  non-homogeneous  processes.  The 
amplitude  of  the  non-homogeneous  components  vary  with  reference  location  and 
therefore  lack  a  comprehensive  characterization.  A  more  comprehensive  characterization 
is  desirable  and  will  be  pursued  in  the  following  chapter. 


*  It  is  noted  ^  that  the  experimental  development  of  an  indicator  function  signal  requires 
the  experimenter  to  set  a  "threshold"  amplitude  level  of  sensor  output  which  identifies 
turbulent  activity.  Errors  in  setting  this  threshold  level  would  affect  the  accuracy  of  spot 
leading  and  trailing  edge  velocities  determined  from  the  indicator  function  signal  only. 
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CHAPTER  3 


TWO-WAVENUMBER-FREQUENCY  SPECTRUM  OF 
TRANSITION  FLOW  PRESSURE  FIELD 


Chapter  2  provides  identification  and  interpretation  of  transition  flow  pressure 
field  features  using  the  space-varying  wavenumber-frequency  spectrum.  This  spectral 
format  allows  immediate  classification  of  dominant  pressure  field  components.  However, 
because  the  non-homogeneous  pressure  field  spectral  amplitude  varies  with  reference 
position,  a  composite  description  of  the  field  cannot  be  obtained.  A  solution  to  this 
problem  is  potentially  available  from  a  spectral  format  defined  by  Strawderman  and 
denoted  the  "Two  Wavevector-Frequency  Spectrum".  This  spectrum,  written  in  terms  of 
the  wavenumber-frequency  spectrum  of  interest  to  this  analysis,  is  defined  as 


(9) 


where  !2^p(77,,£i,T)  is  the  space-time  pressure  field  correlation  and  l^\  and  ^'i  are  the 
Fourier  conjugate  variables  of  the  absolute  spatial  vector  components  rj,  and  £, 
respectively. 


Example  Two-wavenumber  Spectrum 
In  order  to  develop  a  basis  for  the  interpretation  of  the  transition  flow 
pressure,  two  wavenumber-frequency  spectrum,  it  is  useful  to  construct  an  example  of  a 
non-homogeneous  pressure  field.  The  particular  example  selected  is  for  an  exponentially 
decaying  wave  and  the  development  parallels  that  of  Strawderman 
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Let  the  selected  pressure  field  be  defined  as 

p(77,,r)  =  sin 77,  -£Uor  +  0)  (10) 

where  r|^  is  absolute  spatial  position  as  before,  t  is  time,  a  is  a  positive  constant,  and 
COq  are  selected  wavenumber  and  radian  frequency,  and  0  is  a  random  variable  uniformly 
distributed  in  (-k  ,  n).  For  a  =  5,  the  spatial  decay  of  this  field  is  shown  in  Figure  8. 

Note  that  the  amplitude  is  substantially  reduced  over  a  length  consistent  with  that  of  the 


rii  (m) 


Figure  8.  Spatial  Decay  Characteristics  of  Example  Harmonic  Wave 


transition  zone.  The  space-time  correlation  of  this  pressure  field  is 
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-co^t+e) 


(11) 


""■•sin(jto(T7i  +e^)-cOo{t  +  T)+  e)d9 

Performing  the  integration  provides  the  space-time  correlation  of  this  non-homogeneous 
pressure  field 


=  -a)oT)  (12) 

Using  (9),  the  true  two  wavevector-frequency  spectrum  is  given  by 
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From  analysis  of  this  functional  form  Strawderman  notes  that  this  two  wavevector- 
frequency  spectrum  "has  relative  maxima  along  the  lines  k  =  kQ  and  and  an 

absolute  maximum  at  the  intersection  of  these  lines:  that  is  ,  at  ()i,/:)=(0,/:o)"- 

To  aid  visualization  of  this  characteristic  shape  it  is  useful  to  graphically 
present  the  two  wavevector-frequency  spectrum  of  (13).  Also,  to  retain  focus  on 
transition  flow  pressure  fields,  parameters  are  selected  consistent  with  those  of  the  space- 
varying  wavenumber  spectrum  computed  in  Chapter  2.  Therefore  select  cOq  =  306.8 
(rad/s) ,  ito  =  37.1  (rad/m)  and  a  =  5  (m’l).  A  surface  plot  of  (13),  using  these 
parameters,  is  shown  in  Figure  9. 
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Figure  9.  Surface  Plot  of  Exponentially  Decaying  Single  Harmonic  Wave  True  Two- 

Wavenumber  Spectrum 

A  contour  plot  of  these  same  values  is  shown  in  Figure  10.  Note  the  shape  of 
the  contours  in  the  region  of  the  peak  value  and  also  that  the  peak  value  lies  at  the 
intersection  of  the  lines  =  38.13  (rad/m)  and  k-\i=kQ. 

17 

Kennedy  and  Slrawderman  computed  numerical  estimates  of  the  two 
wavenumber-frequency  spectrum  for  an  exponentially  damped  wave.  The  contour  plot  of 
their  estimated  two  wavenumber  spectrum  is  shown  in  Figure  1 1.  Again  note  the  shape 
of  the  contours  in  the  region  of  the  peak  value.  Kennedy  and  Strawderman  call  this  shape 
the  "footprint"  of  an  exponentially  damped  wave. 
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Figure  10.  Contour  Plot  of  Exponentially  Decaying  Single  Harmonic  Y/ave 

Two  Wavenumber  Spectrum 

A  study  of  the  Kennedy  and  Strawderman  work  develops  intuition  for  the 
nature  of  the  two  wavenumber  spectrum.  This  allows  computation  of  an  estimate  for  the 
two  wavenumber  spectrum  of  (13)  by  using  the  finite  Fourier  transforms  implemented  in 
this  research.  This  is  done  for  both  a  uniform  spatial  weighting  function  and  the 
previously  defined  Taylor  window  so  that  windowing  effects  can  be  observed. 
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A 

where  Spp[}l^,k^,co)  identifies  the  two-wavenumber  spectral  estimate,  and  {rAr|^) 

is  an  additional  spatial  weighting  function.  All  other  variables  are  as  previously  defmed. 
For  the  first  calculation  the  spatial  weights  are  set  to  unity.  A  two-wavenumber  spectral 
estimate  surface  plot  of  the  damped  harmonic  wave  is  shown  in  Figure  12. 

A  contour  plot  of  the  same  estimate  is  shown  in  Figure  13.  Note  that  the 
expected  "footprint"  is  indeed  obtained. 


Figure  12.  Surface  Plot  of  Estimated  Two-Wavenumber  Spectrum  of  Example 
Exponentially  Decaying  Single  Harmonic  Wave  Using  Unit  Spatial  Weighting 
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Figure  13.  Contour  Plot  of  Estimated  Two-Wavenumber  Spectrum  of  Example 
Exponentially  Decaying  Single  Harmonic  Wave  Using  Unit  Spatial  Weighting 

Since,  in  general,  weighting  functions  are  required  to  compute  numerical 
spectral  estimates  it  is  appropriate  to  now  apply  them  to  this  idealized  waveform.  This  is 
implemented  by  replacing  the  previously  applied  uniform  weights  with  Taylor  weights 
and  repeating  the  calculation.  A  surface  plot  of  the  two-wavenumber  spectral  estimate 
using  Taylor  weights  is  shown  in  Figure  14. 
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Figure  14.  Surface  Plot  of  Estimated  Two-Wavenumber  Spectrum  of  Example 
Exponentially  Decaying  Single  Harmonic  Wave  Using  Taylor  Spatial  Weighting 

A  contour  plot  of  the  same  estimate  is  shown  in  Figure  15.  Note  that  the  effect  of  the 
weighting  functions  is  to  broaden  the  response  in  the  region  of  the  peak  and  reduce 
amplitudes  in  regions  removed  from  the  peak.  This  is  entirely  consistent  with  the 
previous  discussion  regarding  the  trade-off  between  main  lobe  width  and  sidelobe  level. 
For  this  single  component  example,  the  selected  window  appears  to  be  poor  choice. 
However,  it  must  be  remembered  that  the  window  was  selected  to  reduce  sidelobe  level 
leakage  so  that  low  level  components  might  be  identified.  The  utility  of  this  window  was 
demonstrated  in  Figure  6,  where  the  low  level  trailing  edge  component  velocities  were 
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observed.  It  is  likely  that  these  low  level  components  would  not  have  been  visible  if  a 
window  with  less  sidelobe  rejection  was  used.  It  is  believed  that  this  window  remains  an 
appropriate  choice  for  the  spectral  estimates  of  interest  in  this  research. 

This  example  provides  a  basis  for  the  interpretation  of  the  two- wavenumber 
spectrum  of  a  single,  exponentially  decaying  wave.  It  is  observed  that  this  two- 
wavenumber  spectrum  has  a  characteristic  shape  in  the  ^i,k  spectrum  that  contains  an 
absolute  peak  at  (p.,k)  =  (0,  kg)  and  relative  maxima  along  the  line  k-p,  =  kg. 


k  (rad/m) 

Figure  15.  Contour  Plot  of  Estimated  Two-Wavenumber  Spectrum  of  Example 
Exponentially  Decaying  Single  Harmonic  Wave  Using  Taylor  Spatial  Weighting 
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Two-Wavenumber  Transitional  Flow  Pressure  Field  Spectnim 

An  estimate  of  the  streamwise  transitional  flow  pressure  field  two  wavevector-frequency 
spectrum  is  implemented  by  substituting  into  (14)  yielding 

^  R  N  M 

Spp{^„k,,co)^^  X  Z  ^  {rAr|^)w^^  {nAe^)w^{mAT)R,[rA^J^,nAe^,Q,mAt) 

r=l  n=l  m=l  (25J 

-iiu.rAn,  +  k.nAe.  -  comAx] .  .  , 

e  ^  'Ari^Ae^Ax 

where  R,[r|^,£^,0,x)  is  the  previously  defined  space-time  correlation  of  the  intermittency 
indicator  function,  and  /t,  and  k^  are  now  the  Fourier  conjugate  variables  of  r|^  and  £^ 
respectively. 

The  surface  plot  of  the  50  Hz  two-wavenumber  spectrum  is  shown  in 
Figure  16.  The  contour  plot  of  the  same  spectrum  is  shown  in  Figure  17.  These  spectra 
were  calculated  using  the  Taylor  temporal  and  spatial  weights.  Three  observations  may 
be  made:  1)  the  absolute  spectral  peak  occurs  at  p.  =0  and  the  convective  wavenumber  ( k^ 
=  37.1),  2)  the  line  k-p=A:g,  is  perceptible  only  in  the  region  of  the  peak,  3)  relative 
maxima  occur  at  all  k-wavenumbers  at  |i  =0.  The  location  of  the  absolute  peak  is  as 
expected  and  needs  no  further  discussion.  However,  the  other  two  observations  warrant 
further  consideration. 

Recall,  from  the  exponentially  decaying  wave  example,  that  the  effect  of  the 
application  of  spatial  weighting  was  to  obscure  the  visualization  of  the  relative  maxima 
line  (k-p=  k^ ).  This  is  likely  the  result  of  the  wider  main  lobe  and  lower  sidelobes  of  the 

spatial  window.  Also,  recall  that  the  example  was  for  a  single  wave  at  an  arbitrary 
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Figure  16.  Surface  Plot  of  Estimated  Streamwise  Transitional  Flow  Pressure  Field 

Two-Wavenumber  Spectrum  at  50  Hz 

wavenumber  .  It  is  to  be  expected  that  another  independent  wave,  at  a  different  k- 
wavenumber,  would  create  a  similar  shape  but  the  footprint  would  be  displaced 
somewhere  else  consistent  with  the  selected  wavenumber.  Now  assume  that  the 
transitional  flow  pressure  field  is  the  superposition  of  a  continuous  distribution  of  such 
components.  With  this  assumption  the  wavenumber  spectra  at  |i  =0  would  represent  the 

propagating  wavenumber  spectrum  of  the  transition  flow  pressure  field.  Furthermore,  the 
spectral  contributions  of  any  given  component  along  the  line  k-|i=^o  would  merge  with 
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Figure  17.  Contour  Plot  of  Estimated  Streamwise  Transition  Flow  Pressure  Field 

Two-Wavenumber  Spectrum  at  50  Hz 

those  of  adjacent  spectral  components,  destroying  the  characteristic  "footprint"  observed 
in  the  example.  With  these  considerations  it  is  believed  that  Figures  16  and  17  are  valid 
representations  of  the  two-wavenumber  spectrum  of  the  transitional  flow  pressure  field. 

To  complete  this  presentation  the  two-wavenumber  spectrum  at  100  Hz  is 
calculated.  The  surface  plot  of  this  spectrum  is  shown  in  Figure  18. 
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Figure  18.  Surface  Plot  of  Estimated  Streamwise  Transition  Flow  Pressure  Field 

Two-Wavenumber  Spectrum  at  100  Hz 

The  contour  plot  is  shown  in  Figure  19.  Note  that  the  absolute  peak  value  is  again  at  the 
convective  wavenumber.  Also  note  that  a  lower  but  significant  peak  occurs  at  a  k 
wavenumber  of  approximately  -116.  This  is  consistent  with  the  turbulent  spot 
generation  process  noted  in  the  space-varying  spectrum.  It  is  apparent  then  that  the  two- 
wavenumber  spectrum  provides  the  same  information  available  from  the  multiple  space- 
varying  wavenumber  spectra,  but  in  a  single  concise  format. 
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Two-Wavenumber  Spectrum  at  100  Hz 


Summary  of  the  Two-Wavenumber  Spectra 
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A  two-wavenumber  spectrum  has  been  defined  that  provides  a  concise 
description  of  non-homogeneous  fields.  The  characteristics  of  this  two- wavenumber 
spectrum  have  been  identified  using  a  non-homogeneous  pressure  field  example 
containing  parameters  pertinent  to  a  subsonic  boundary  layer  transition  zone.  Transition 
zone  flow  pressure  field  estimates  for  the  two- wavenumber  spectrum  have  been 
computed  at  two  frequencies  of  interest.  These  two-wavenumber  spectral  estimates 
provide  a  full  spectral,  concise  form  of  the  information  contained  in  the  associated 
multiple  space-varying  wavenumber-frequency  spectra.  The  results  of  this  research 
demonstrate  that  the  two-wavenumber  spectrum  provides  a  concise,  full  spectral 
representation  of  the  streamwise  transitional  flow  pressure  field. 


CHAPTER  4 
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CONCLUSIONS  AND  RECOMMENDATIONS 


This  research  has  implemented  a  computer  model  to  produce  unprecedented 
images  of  the  non-homogeneous,  transitional  flow  pressure  field  wavenumber-frequency 
spectrum,  in  two  formats.  The  first  format,  called  the  space-varying  wavenumber- 
frequency  spectrum,  accommodates  non-homogeneous  fields  by  presenting  the  familiar 
wavenumber-frequency  spectrum  in  terms  of  an  absolute  spatial  reference  position 
parameter.  Because  the  field  is  non-homogeneous,  multiple  wavenumber  spectra  must  be 
computed  to  display  the  field  characteristics.  This  research  provides  these  multiple 
wavenumber  spectra  at  two  frequencies  of  interest.  Analysis  of  these  space-varying 
wavenumber-frequency  spectral  estimates  suggest  that  more  than  one  convective 
component  exists  in  the  transition  zone  indicator  function. 

There  is  a  consistently  prominent  convection  wavenumber  where  JJ^ 

is  the  average  convection  speed  of  the  turbulent  spots.  And,  there  also  appears  a  higher 
wavenumber  component  corresponding  to  a  much  lower  convection  speed  than  the 
average  convection  speed.  The  speeds  computed  from  the  spectral  information  suggest 
that  these  lower  convection  speeds  are  associated  with  the  much  slower  trailing  edge 
velocity  of  the  turbulent  spots  within  the  U'ansition  zone.  This  identification  would  not 
have  been  possible  from  a  study  of  the  space-time  correlation  functions  alone.  This 
demonstrates  the  power  of  the  wavenumber  spectral  representation  in  uncovering 
important  physical  features  of  complex  non-homogeneous  fields. 


37 

The  second  wavenumber-frequency  spectral  format  is  newly  defined  and  called 
the  two-wavenumber-frequency  spectrum.  This  format  is  an  extension  of  the  space- 
varying  wavenumber-frequency  spectrum  and  provides  a  full  spectral  representation  of  a 
non-homogeneous  field.  This  format  has  the  distinct  advantage  of  providing  a  concise 
presentation  of  the  wavenumber  spectral  properties  of  the  field,  at  a  given  frequency,  in  a 
single  image.  Associated  with  the  advantages  of  conciseness  is  the  need  to  develop 
knowledge  and  intuition  with  the  new  and  somewhat  unfamiliar  jj,  wavenumber  domain. 

It  is  recommended  that  this  knowledge  be  pursued  and  that  this  spectral  technique 
be  utilized  in  future  experimental  efforts  involving  non-homogeneous  random  fields^  It 
would  be  particularly  valuable  in  transition  zone  wall  pressure  fluctuation  studies  to 
provide  direct  measurement  of  the  pressure  field.  With  the  entire  pressure  signal 
available,  it  is  believed  that  the  two-wavenumber  spectral  representation  would  yield 
accurate  estimates  of  the  leading  and  trailing  edge  spot  convection  velocities  as  well  as 
the  average  spot  convection  velocity.  Remember  that  the  indicator  function  used  here 
was  for  natural  turbulent  spots  that  are  created  at  random  in  both  space  and  time.  The 
work  of  Ref  [15],  which  identified  leading  and  trailing  edge  celerities,  used  single  spots 
created  deterministically.  Here,  we  have  uncovered  almost  the  same  information  for 
randomly  occurring  spots  in  a  naturally  occurring  transition  zone. 
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